International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-3, Issue-3 March- 2017] 


New Finite Element for modeling connections in 2D small frames 

- Static Analysis 

Chadi Azoury 1 , Assad Kallassy 2 , Ibrahim Moukarzel 3 

Mechanical Engineering Department, Lebanese University, Lebanon 
2 Civil Engineering Department, Lebanese University, Lebanon 
3 Faculty of Engineering, Lebanese University, Lebanon 


Abstract — The paper presents a proposal to create a new one -dimensional (1-D) finite element that can substitute for the 
three-dimensional (3-D) finite model of connection elements in static analysis. The study is conducted on a simple small 
frame. First, a 3-D model is constructed, followed by the construction of an equivalent two-dimensional ( 2-D ) one. At the 
final stage, the results of the static analysis of many proposed 1-D element models are compared with the results of the 2-D 
model. With a special treatment of the connection elements, the results of the analysis of the new 1-D finite element models 
agree well with the results of the 2-D model. 
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Notations 

E Young’s elastic modulus 

v Poisson coefficient 

p density 

I moment of inertia 

E I flexural beam stiffness = E x I 

t frame thickness 

K(nxn) stiffness matrix of order n 

L (nxm) rectangular matrix with n rows and m columns 
EB Euler-Bernoulli 

sh shear deformation 

FE finite element 

DOF degree-of-freedom 


I. Introduction 

To perform a static analysis of a frame, 3-D finite elements are usually used in the meshing process. A large number of 
elements are required, which in turn requires a large number of DOF and consequently, a great deal of work. The main 
objective consists of developing a 1-D model made up of simple 1-D elements to replace the 3-D model made up of 
tetrahedral prisms or brick elements. To perform a static analysis, the calculation of the stiffness matrix is required. 

In order to find the best 1-D finite meshing element, 3-D meshing is used, and 2-D and 1-D models are consecutively 
constructed. 

All calculations are done by using ANSYS Version 16.2 and MATLAB Version 7. A simple frame is used in the 
calculations. 


II. Test structure 

Our test structure will be of this form: 

0.075 m 


K >1 



Figure 1: Test structure 
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The beam section of the test structure is rectangular (0.02x0.02 m 2 ). The mechanical characteristics are: Poisson coefficient 
v = 0.31, density p = 78000 N/m 3 and elastic modulus E = 184 GPa. The test structure is clamped at its base and consists of 5 
structural elements shown in Figure 2. 



Figure 2: Test structure constructed in ANSYS® R16.2 


Different types of elements are used in the meshing processes. In static analysis, the comparison between different models is 
based on the deformation in the XY plan of the structure. These models are compared in 2 load cases illustrated in 

Figure 3: 

• Case 1 : a body load of 78000 N/m 3 along the vertical axis of the horizontal beam 

• Case 2: a body load of 78000 N/m 3 along the horizontal axis of the left column 




No load 
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Body load along X + 
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Central line \ 

Load case 2 
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Figure 3: Test structure with load cases 1 and 2 
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III. 3D AND 2D MODELING 


The 3D model of the test structure is meshed using H8 volumetric finite elements with 8 nodes. Each column is modeled 
using 8x8x16 = 1024 elements and each corner includes 8x8x8 = 512 elements. The horizontal beam is modeled using 
8x8x16 = 1024 elements. In total, it is equal to 4096 elements with 5265 nodes (see Figure 4). 



Figure 4: 3D test structure modelled using H8 Finite Elements 


Figures 5 and 6 show the total displacement contours of the frame under load case 1 and 2. 



Figure 5: Three-dimensional deformed test structure under load case 1 



Figure 6: Three-dimensional deformed test structure under load case 2 
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The test structure is then modeled in two dimensions. A Q4 plane element with a thickness t = 0.02 m and having 4 nodes is 
used. Each column is modeled using 8x16 = 128 elements. Each corner includes 8x8 = 64 elements. The horizontal beam is 
modeled using 8x16 = 128 elements. The total is 512 elements with 585 nodes (see Figure 9a). 



Figure 7: 2D Test structure modeled using Q4 plane Finite Elements 
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Figure 8: Two-dimensional deformed test structure along load case 1 



Figure 9: Two-dimensional deformed test structure along load case 2 
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Figures 8 and 9 have similar maps as Figures 5 and 6. The different models are compared using the displacement of 65 points 
on the central line of the longitudinal symmetrical plan of the structure. 



Figure 10 : Longitudinal symmetrical plan of the structure showing the 65 central line 

points 

In order to establish a comparison between both FE models and using both load cases, UX and UY displacements of nodes 
on the central line, along the X and Y axis respectively, are plotted on the same graph (see Figure 11-10). 



Nodes 

Figure 11 : UX (m) central line displacement along axe X (load case 1) 



Nodes 

Figure 12 : UY (m) central line displacement along axe Y (load case 1) 
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Nodes 

Figure 13 : UX (m) central line displacement along axe X (load case 2) 



Nodes 

Figure 14 : UY (m) central line displacement along axe Y (load case 2) 

Regarding the comparison criteria of the results produced by different models, the following error formula is used: 

\U R (i)-U c (i)\ 


E mm - max 


ma x(U R (;')) - min^ (;)) 


U R (i) = reference value of a displacement at point i 

U c (i) = calculated value of the same displacement and at the same point i of a new model 
The maximum errors for both models H8 and Q4, for both load cases, are listed in Table 1 below: 


Table 1 

Comparison between H8 and Q4 models 


Load case 

Displacement 

Utf(m) 

[H8] 

Uc(m) 

[Q4] 

Emax 

Node number (Emax) 

1 

UX 

-2.6673E-10 

-2.7536E-10 

1.55% 

15 and 51 

1 

UY 

-1.0941E-09 

-1.1149E-09 

1.90% 

33 

2 

UX 

2.4111E-09 

2.4842E-09 

3.01% 

16 

2 

UY 

2.6802E-10 

2.7866E-10 

2.92% 

39 


The plane model has produced comparable results to those of the 3-D model. Thus, the reference model is the 2-D model 
which leads to a reduction in the DOF amount. It will be used to develop the best 1-D element for connection modeling. 
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IV. One-dimensional classical modeling 

In this section, 1-D elements are used for the test structure. Since the structure is not the slender type, we know that Euler- 
Bernoulli beams are not adequate. Therefore, we will use Timoshenko beams with 2 nodes per element and 3 DOF per node 
(Timoshenko model). 
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Figure 15 : Two-node beams 

The transformation from the Q4 model to the 1-D model is shown in Figures 12 and 13. The volumetric charge is 
transformed into linear charge. 


78 KN/m 3 





Figure 16 : Transformation of the Q4 model into a 1-D model (load case 1) 
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Figure 17 : Transformation of the Q4 model into a 1-D model (load case 2) 


The left and right columns are modeled using 16 elements each. The horizontal beam is also modeled using 16 elements. 
Both left and right corners are modeled using 4 horizontal linear elements and 4 vertical linear elements, i.e. a total of 64 
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beam elements. The following figures show the deformed test structure in the 2 load cases as applicable to both models: the 
reference model Q4 and the Timoshenko model. 



Nodes 

Figure 18 : UX (m) displacement of central line nodes along X (load case 1) 



Figure 19 : UY (m) displacement of central line nodes along Y (load case 1) 



Nodes 


Figure 20 : UX (m) displacement of central line nodes X (load case 2) 
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Figure 21 : UY (m) displacement of central line nodes Y (load case 2) 


None of the 1-D models produces proximate results in comparison with those given by the Q4 model. Improved models 
using Timoshenko beams are shown below. 

V. First improved model: CQ corner quad 

A new CQ element is suggested for the corners, based on the utilization of Q4 surface elements since the Q4 modeling and 
the H8 modeling produce close results. Non-corner elements are still modeled as Timoshenko beams. One way to connect 
two different meshing types is explained by Dohrmann and al, [4] and [5], and Quiroz and al [9]. Kattner and al [8] explain 
how modeling joints between frame columns and beams is performed. 

The interface between the new CQ element on the one hand and the horizontal and vertical beams on the other hand, is 
considered rigid in our case. The corner element of 0.02 m x 0.02 m and 0.02 m in thickness is meshed using Q4 elements of 
0.02/8 m x 0.02/8 m which produces 64 elements and 81 nodes in total, hence 81x2 DOF: 17 x 2 interface DOF and 64 x 2 
non-interface DOF (see next Figure). 


D 



Figure 22 : Left corner with rigid links 

The elementary stiffness matrix kg(8,8) of the Q4 element is mentioned in Cook [2]. By assembling the 64 Q4 elements, the 
global stiffness matrix of the corner is K(162,162). 


Using static condensation, the global stiffness matrix is reduced to K c (34,34) after elimination of the 128 non-interface DOF 
[3]. 
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To find matrix Kc(34,34), we use the formula for static condensation. In Matlab code: 

krr=K(l:34, 1:34) ; 
krc=K(l:34, 35:162) / 
kcr=K (35:162,1:34) / 
kcc=K (35:162,35:162) ; 

Kc=krr-krc*inv ( kcc) *kcr; 

The AB interface is supposed to be rigid. Then, 3 DOF (q b q 2 , q 3 ) are enough to determine the DOF of the AB interface: 

qi: horizontal displacement of the midpoint of segment AB 
q 2 : vertical displacement of the midpoint of segment AB 
q 3 : rotation of segment AB around its midpoint 

The same applies to the BC interface where the DOF (q 4 , q 5 , q 6 ) of the midpoint in BC determine the DOF of the rigid 
interface BC. 

The reduced stiffness matrix K r (6,6) is given by: K r = L T K C L, where L(34x6) is the matrix giving the values of the DOF 
relative to interfaces AB and BC, from degrees qi, q 2 , q 3 , q 4 , q 5 , and q 6 . 

L=[ 

1 0 0 0 0 0 

0 1 -4 *a 0 0 0 

1 0 0 0 0 0 

0 1 -3*a 0 0 0 

1 0 0 0 0 0 

0 1 -2 *a 0 0 0 

1 0 0 0 0 0 

0 1 -l*a 0 0 0 

1 0 0 0 0 0 

0 1 0 0 0 0 

1 0 0 0 0 0 

0 1 l*a 0 0 0 

1 0 0 0 0 0 

0 1 2 *a 0 0 0 

1 0 0 0 0 0 

0 1 3*a 0 0 0 

0 0 0 1 0 4*b 

0 1 4 *a 0 0 0 

0 0 0 1 0 3*b 

0 0 0 0 1 0 

0 0 0 1 0 2 *b 

0 0 0 0 1 0 

0 0 0 1 0 l*b 

0 0 0 0 1 0 

0 0 0 1 0 0 

0 0 0 0 1 0 

00 0 1 0 -l*b 

0 0 0 0 1 0 

0 0 0 1 0 -2 *b 

0 0 0 0 1 0 

00 0 1 0 -3*b 

0 0 0 0 1 0 

0 0 0 1 0 -4 *b 

0 0 0 0 1 0 ] / 

The stiffness matrix K r (6x6) is a function of: 

• The length a of side AB 

• The length b of side BC 
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• The elastic modulus E of the material in part ABCD 

• The Poisson coefficient v 

• The thickness of the frame t 

The stiffness matrix of the left corner beam of the test structure is given numerically by: 

kief t— [ 

4 . 6660e+00 9 1.5022e+009 -5.2058e+006 -4.6660e+009 -1.5022e+009 -2.6432e+007 

1 . 5022e+00 9 4.6660e+009 2.6432e+007 -1.5022e+009 -4.6660e+009 5.2058e+006 

-5 . 2058e+006 2.6432e+007 3.7398e+005 5.2058e+006 -2.6432e+007 -5.7599e+004 

-4 . 6660e+00 9 -1.5022e+009 5.2058e+006 4.6660e+009 1.5022e+009 2.6432e+007 

-1 . 5022e+00 9 -4.6660e+009 -2.6432e+007 1.5022e+009 4.6660e+009 -5.2058e+006 

-2 . 6432e+007 5.2058e+006 -5.7599e+004 2.6432e+007 -5.2058e+006 3.7398e+005] 

With a simple transformation, we can obtain the stiffness matrix of the right corner beam: 

kright= [ 

4 . 6660e+00 9 -1.5022e+009 -5.2058e+006 -4.6660e+009 1.5022e+009 -2.6432e+007 

-1 . 5022e+00 9 4.6660e+009 -2.6432e+007 1.5022e+009 -4.6660e+009 -5.2058e+006 

-5 . 2058e+006 -2.6432e+007 3.7398e+005 5.2058e+006 2.6432e+007 -5.7599e+004 

-4 . 6660e+00 9 1.5022e+009 5.2058e+006 4.6660e+009 -1.5022e+009 2.6432e+007 

1 . 5022e+00 9 -4.6660e+009 2.6432e+007 -1.5022e+009 4.6660e+009 5.2058e+006 

-2 . 6432e+007 -5.2058e+006 -5.7599e+004 2.6432e+007 5.2058e+006 3.7398e+005] 

Finally, the CQ model of the constructed frame is as follows: 



Figure 23 : Description of the CQ model 


MN, OP and QR are Timoshenko beams 

NO is a CQ element with a K rleft (6x6) stiffness matrix 

PQ is a CQ element with a K rright (6x6) stiffness matrix 

Figures 24 to 27 show static results under load cases 1 and 2. 
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Figure 25 : UY (m) displacement of the central line nodes along Y (load case 1) 



Figure 26 : UX (m) displacement of the central line nodes along X (load case 2) 
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Figure 27 : UY (m) displacement of the central line nodes along Y (load case 2) 

The graphs show that the results of the CQ model are very close to those of the Q4 model. Table 2 shows the results of both 
models for both load cases. The last column gives the node number corresponding to the highest error value E max . 


Table 2 

Comparison of both models and all load cases 


Load 

case 

Displacement 

U*(m) 

[Q4] 

Model 

U c (m) 

Emax 

Node number 

( Emax) 

1 

UX 

2.8613E-10 

Timoshenko 

1.9284E-10 

16.28% 

14 or 52 

1 

UY 

-1 .1 149E-09 

Timoshenko 

-1.5148E-09 

35.87% 

33 

2 

UX 

2.3397E-09 

Timoshenko 

2.7944E-09 

18.19% 

36 or 37 

2 

UY 

2.7866E-10 

Timoshenko 

1.645E-10 

30.38% 

39 

1 

UX 

-2.7536E-10 

CQ 

-2.6053E-10 

2.59% 

15 or 51 

1 

UY 

-3.1495E-10 

CQ 

-3.3374E-10 

1.69% 

17 or 49 

2 

UX 

2.3947E-09 

CQ 

2.4623E-09 

2.70% 

28 

2 

UY 

9.7806E-1 1 

CQ 

8.2391 E-11 

4.10% 

25 


It can be said from the previous table, that the CQ model produces the best approximation to the Q4 model. 

VI. Second improved model: Rigid interface everywhere (RQ) 


In this paragraph, we propose a linear finite element to replace Timoshenko beams in the vertical and horizontal beams. This 
new element is constructed using the same procedure for the corner elements as in the previous paragraph. 

Consider a strip in the beam that is between two consecutive nodes of the midline (see next figure). 




a 



Figure 28: RQ element 


Page 1 114 


International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-3, Issue-3 March- 2017] 


The total length of this element is 0.02 m (from bottom to top). The thickness (along z direction) is always 0.02 m. a = b = 
0.045/16 = 0.0028125 m. 

The element is discretized in 8 Q4 elements along its length. Therefore, we have a total of 9 nodes and 18 DOF for each 
interface. 

We will assume that plane sections of beams in bending (Timoshenko hypothesis), displacements of nodes that are in the 
same cut-section are displacements of a rigid body, given by 3 DOFs: 

• UX = Displacement along the X axis of the center of the section 

• UY = Displacement along the Y axis of the center of the section 

• RZ = Rotation around the Z axis 

Therefore, the displacements of all the nodes of this element are given by 6 DOFs; 3 DOFs for the left cut-section and 3 
DOFs for the right cut-section. This will give us the elementary stiffness matrix ke(6,6). In order to determine ke(6,6), we 
know that the stiffness matrix for one Q4 element is given by: 

The global matrix is K(36,36). 

The matrix of the rigid interface is given by: 

R=[ 

1 0 4 *b 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 4 *b 

0 0 0 0 1 0 

1 0 3*b 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 3*b 

0 0 0 0 1 0 

1 0 2 *b 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 2 *b 

0 0 0 0 1 0 

1 0 l*b 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 l*b 

0 0 0 0 1 0 

1 0 0 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 0 

0 0 0 0 1 0 

1 0 -l*b 0 0 0 

0 1 0 0 0 0 

00 0 1 0 -l*b 

0 0 0 0 1 0 

1 0 -2 *b 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 -2 *b 

0 0 0 0 1 0 

1 0 -3*b 0 0 0 

0 1 0 0 0 0 

00 0 1 0 -3*b 

0 0 0 0 1 0 

1 0 -4 *b 0 0 0 

0 1 0 0 0 0 

0 0 0 1 0 -4 *b 

0 0 0 0 1 0 ] 

Therefore, K = L T xKcxL. (Kc is the condensed matrix). We obtain: 


Page 1 115 




International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-3, Issue-3 March- 2017] 


k=E*t/ ( v A 2 -1 ) *b/a* [ 

- 8 , 0 , 0 , 8 , 0 , 0 

0, 4 *v-4 , 2*a*(v-l), 0, -4*v+4, 2*a*(v-l) 

0 , 2*a* (v-1) , -128/3*b A 2-4/3*a A 2+4/3*a A 2*v, 0,-2*a* (v-1) , 128/3*b A 2-2/3*a A 2+2/3*a A 2*v 

8 , 0 , 0 ,- 8 , 0 , 0 

0, -4 *v+4 , -2*a*(v-l), 0, 4*v-4, -2*a*(v-l) 

0 , 2*a* (v-1) , 128/3*b A 2-2/3*a A 2+2/3*a A 2*v, 0,-2*a* (v-1) , -128/3*b A 2-4/3*a A 2+4/3*a A 2*v] 

Our new modeling is shown in the next figure: 



Figure 29: RQ model for the horizontal and vertical beams while CQ element at the corners 

HO and O’H’ are beam elements formed from 8 Q4 elements. OO’ is a corner element. 

The following graphs show the displacements along X and Y, for the two load cases for all the models that we have 
considered so far (Q4, Timoshenko, CQ, and RQ): 



Figure 30 : UX (m) displacement of the central line nodes along X (load case 1) 


Page I 116 


International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-3, Issue-3 March- 2017] 



Figure 31 : UY (m) displacement of the central line nodes along Y (load case 1) 



Figure 32: UX (m) displacement of the central line nodes along X (load case 2) 



Figure 33 : UY (m) displacement of the central line nodes along Y (load case 2) 
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VII. Fourth improved model (RQw) 

The graphs of the previous paragraph show that the RQ model is a little stiffer than the CQ model. In order to give it a certain 
flexibility, we will add two DOFs (warping) for every node of the previous paragraph: 

A new element is constructed from the RQ model that allows a plane section of a beam to deform a cubic or quadratic 
deformation after solicitation. Therefore, two DOFs will be added on each node: 

A warping (in terms of y 2 ) that allows a quadratic displacement of the plane section and another warping term (in terms of y 3 ) 
that allows a cubic displacement of the plane section. The following figure illustrates the 5 DOFs on each node for the new 
element. 



Figure 34: The 5 DOFs for the RQw model 

For this element with 5 DOFs per node we expect a stiffness matrix k(10,10). The terms of the 4 th and 5 th degree of the 
stiffness degree are calculated as follows: 

We find the global matrix K(36,36), then matrix L of the rigid link: 


L=[ 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 


4 *b 
0 
0 
0 

3*b 

0 

0 

0 

2*b 

0 

0 

0 

l*b 

0 

0 

0 

0 

0 

0 

0 

-l*b 

0 

0 

0 

-2*b 

0 

0 

0 

-3*b 

0 

0 

0 

-4 *b 
0 
0 
0 


-1 6*b A 2 
0 
0 
0 

-9*b A 2 

0 

0 

0 

-4 *b A 2 
0 
0 
0 

-l*b A 2 

0 

0 

0 

0 

0 

0 

0 

-l*b A 2 

0 

0 

0 

-4 *b A 2 
0 
0 
0 

- 9*b A 2 
0 
0 
0 

-16*b^ 

0 

0 

0 


64 
0 
0 
0 

2 7 *b^ 
0 
0 
0 

8*b A 2 

0 

0 

0 

l*b A 2 

0 

0 

0 

0 

0 

0 

0 

-l*b A 

0 


b A 3 


-8*b A 3 


0 
0 

-6 
0 
0 
0 

-2 7 *b A 3 0 
0 
0 
0 

-64 
0 
0 
0 


0 

1 

0 

^3 0 
0 
1 
0 


0 

0 

4 *b 
0 
0 
0 

3*b 

0 

0 

0 

2*b 

0 

0 

0 

l*b 

0 

0 

0 

0 

0 

0 

0 

-l*b 

0 

0 

0 

-2*b 

0 

0 

0 

-3*b 

0 

0 

0 

-4 *b 
0 


0 

0 

-1 6*b A 2 
0 
0 
0 

- 9*b A 2 
0 
0 
0 

-4 *b A 2 
0 
0 
0 

-l*b A 2 

0 

0 

0 

0 

0 

0 

0 

-l*b A 2 

0 

0 

0 

-4 *b A 2 
0 
0 
0 

-9*b A 2 

0 

0 

0 

-1 6*b A 2 
0 


0 

0 

64*b A 3 

0 

0 

0 

2 7 *b A 3 
0 
0 
0 

8*b A 3 

0 

0 

0 

l*b A 3 

0 

0 

0 

0 

0 

0 

0 

-l*b A 3 

0 

0 

0 

-8*b A 3 

0 

0 

0 

-2 7 *b A 3 
0 
0 
0 

-64*b A 3 

0 


] ? 
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We can find: 

k = l t *k*l 


the stiffness matrix k for the left comer with (warp) is given numerically by: 

kleft = 


4 . 6660e+009 

1 . 5022e+009 

-5 . 2058e+006 

-7 . 3236e+004 

-4 . 3907e+002 

-4 . 6660e+009 

-1.50226+009 

-2 . 6432e+007 

2 . 5186e+005 

-1 . 8964e-003 


1 . 5022e+009 

4 . 6660e+009 

2 . 6432e+007 

-2 . 5186e+005 

1 . 8964e+003 

-1 . 5022e+009 

-4 . 6660e+009 

5 . 2058e+006 

7 . 3236e+004 

4 . 3907e+002 


-5 . 2058e+006 

2 . 6432e+007 

3 . 7398e+005 

-1. 6755e+003 

2 . 7583e+001 

5 . 2058e+006 

-2 . 6432e+007 

-5 . 7599e+004 

-1 . 1072e+002 

-4 . 2280e+000 


-7 . 3236e+004 

-2 . 5186e+005 

-1 . 6755e+003 

2 . 8150e+001 

-1.2750e-001 

7 . 3236e+004 

2 . 5186e+005 

-1 . 1072 e+ 002 

-4 . 7061e+000 

8 . 7138e-005 


-4 . 3907e+002 

1.8964e+003 

2 . 7583e+001 

-1 . 2750e-001 

2 . 6159e-003 

4 . 3907e+002 

-1 . 8964e+003 
-4 . 2280e+000 

8 . 7138e-005 

-4 . 0191e-004 


-4 . 6660e+009 

-1 . 5022e+009 

5 . 2058e+006 

7 . 3236e+004 

4 . 3907e+002 

4 . 6660e+009 

1.5022e+009 

2 . 6432e+007 

-2 . 5186e+005 

1 . 8964e+003 


-1 . 5022e+009 

-4 . 6660e+009 

-2 . 6432e+007 

2 . 5186e+005 

-1 . 8964e+003 

1.5022e+009 

4 . 6660e+009 

-5 . 2058e+006 

-7 . 3236e+004 

-4 . 3907e+002 


-2 . 6432e+007 

5 . 2058e+006 

-5 . 7599e+004 

-1 . 1072e+002 

-4 . 2280e+000 

2 . 6432e+007 

-5 . 2058e+006 

3 . 7398e+005 

-1. 6755e+003 

2 . 7583e+001 


2 . 5186e+005 

7 . 3236e+004 
— 1.107 2e+0 02 
-4 . 7061e+000 

8 . 7138e-005 

-2 . 5186e+005 

-7 . 3236e+004 

-1.6755e+003 

2 . 8150e+001 

-1 . 2750e-001 


-1 . 8964e+003 

4 . 3907e+002 

-4 . 2280e+000 

8 . 7138e-005 

-4 . 0191e-004 

1.8964e+003 

-4 . 3907e+002 

2 . 7583e+001 

-1 . 2750e-001 

2 . 6159e-003 


With a simple transformation, we can find the stiffness matrix of the right corner of our structure: 

kright = 


4 . 6660e+009 

-1.50226+009 
-5 . 2058e+006 

-7 . 3236e+004 

-4 . 3907e+002 

-4 . 6660e+009 

1.5022e+009 

-2 . 6432e+007 

2 . 5186e+005 

-1 . 8964e-003 


-1. 5022e+009 
4 . 6660e+009 

-2 . 6432e+007 

2 . 5186e+005 

-1 . 8964e+003 

1 . 5022e+009 

-4 . 6660e+009 

-5 . 2058e+006 

-7 . 3236e+004 

-4 . 3907e+002 


-5 . 2058e+006 

-2 . 6432e+007 

3 . 7398e+005 

-1.6755e+003 

2 . 7583e+001 

5 . 2058e+006 

2 . 6432e+007 

-5 . 7599e+004 

■1 . 1072e+002 

-4 . 2280e+000 


■7 . 3236e+004 

2 . 5186e+005 

■1 . 67 55e+003 

2 . 8150e+001 

■1 . 2750e-001 

7 . 3236e+004 

■2 . 5186e+005 

■1 . 1072e+002 

•4 . 7061e+000 

8 . 7138e-005 


-4 . 3907e+002 

-1 . 8964e+003 

2 . 7583e+001 

*1, 2750e-001 

2 . 6159e-003 

4 . 3907e+002 

1. 8964e+003 

-4 . 2280e+000 

8 . 7138e-005 

-4 . 0191e-004 


•4 . 6660e+009 

1 , 5022e+009 

5 . 2058e+006 

7 . 3236e+004 

4 . 3907e+002 

4. 6660e+009 

■1 . 5022e+009 

2 . 6432e+007 

■2 . 5186e+005 

1.8964e+003 


1 . 5022e+009 

-4 . 6660e+009 

2 . 6432e+007 

-2 . 5186e+005 

1 . 8964e+003 

-1 . 5022e+009 

4 . 6660e+009 

5 . 2058e+006 

7 . 3236e+004 

4 . 3907e+002 


■2 . 6432e+007 

•5 . 2058e+006 

•5 . 7599e+004 

■1 . 1072e+002 

•4 . 2280e+000 

2 . 6432e+007 

5 . 2058e+006 

3 . 7398e+005 

■1 . 67 55e+003 

2 . 7583e+001 


2 . 5186e+005 

-7 . 3236e+004 

-1 . 1072e+002 

-4 . 7061e+000 

8 . 7138e-005 

-2 . 5186e+005 

7 . 3236e+004 

-1. 6755e+003 

2 . 8150e+001 

-1 . 2750e-001 


-1.8964e+003 

-4 . 3907e+002 

-4 . 2280e+000 

8 . 7138e-005 

-4 . 0191e-004 

1 . 8964e+003 

4 . 3907e+002 

2 . 7583e+001 

-1 . 2750e-001 

2 . 6159e-003 


The following figures show the displacements along X et Y, for both the load cases, for all models (Reference model: Q4 
model, Timoshenko model, CQ model, RQ, model and RQw). 



Figure 35: UX (m) displacement of the central line nodes along X (load case 1) 
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Figure 36: UY (m) displacement of the central line nodes along Y (load case 1) 



Figure 37: UX (m) displacement of the central line nodes along X (load case 2) 



Figure 38: UY (m) displacement of the central line nodes along Y (load case 2) 
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The previous graphs show that the RQw is more flexible than RQ, whereas CQ gives best results. 

VIII. Summary 

The following table summarizes the results coming from all models: 


Table 3 

Comparison of all models and all load cases 


Model 

Node 

(where E = E max ) 

Load case 

Displacement 

^max 

par rapport a H8 

^max 

par rapport a Q4 

HR 

IE on El 

1 

I IX 

o% 

1 EE% 

HR 

33 

1 

NY 

o% 

1 90% 

HR 

Ifi 

9 

l IX 

o% 

3 01% 

HR 

39 

9 

NY 

o% 

9 99% 

04 

IE on El 

1 

NX 

1 EE % 

0% 

04 

33 

1 

1 IY 

1 90% 

0% 

04 

Ifi 

9 

IJX 

3 01% 

0% 

04 

39 

9 

l IY 

9 99% 

0% 

Timoshenko 

14 on E9 

1 

l IX 

IE 99% 

Ifi 93% 

Timoshenko 

33 

1 

IJY 

3R 4E% 

3E R7% 

Timoshenko 

3fi ON 37 

9 

l IX 

91 43% 

1 R 19% 

Timoshenko 

3R on 39 

9 

NY 

9R 44% 

30 33% 

no 

IE on El 

1 

NX 

1 E9% 

9 EQ% 

no 

17 nil 49 

1 

LJY 

9 74% 

1 R9% 

no 

9R 

9 

NX 

E E9% 

9 70% 

no 

9E 

9 

NY 

4 4R% 

4 in<>A 

RO 

14 ON E9 

1 

IJX 

fi 30% 

7 fi1% 

RO 

33 

1 

NY 

7 71% 

9 44% 

RO 

Ifi 

9 

IJX 

fi 4R% 

9 19% 

RO 

37 

9 

IJY 

R 79% 

10 97% 

ROw 

14 on E9 

1 

NX 

4 34% 

E 71% 

ROw 

33 

1 

IJY 

4.34% 

E 37% 

ROw 

17 

9 

IJX 

1 fifi% 

4 E9% 

ROw 

36 

2 

L1Y 

6 77% 

7 84% 


It is obvious that the CQ model gives best results when comparing with the Q4 model. 

IX. Conclusion 

Replacing the volumetric modeling of a simple frame by traditional 1-D elements leads to a weak simulation in static 
analysis. It is found out that the problem arises from corner connections. Modeling these connections by rigid links produces 
a stiff structure. Using CQ elements at the corners improves the results in the static analysis. Consequently, the number of 
elements and the number of DOF of the finite element model are dramatically reduced from hundreds of thousands to merely 
hundreds. 
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